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^ ' Abstract 



O 



In this work, we present the analytic treatment of the Fourier Trans- 
form (FT) of the Gravitational Wave (GW) signal from a pulsar including 
spin down corrections in a parametrized model discussed by Brady et. al. 
pp. The formalism lends itself to a development of the FT in terms of well 
known special functions and integrals defining the spin down moments. 



> 

Tj- ; 1 Introduction 

' The detection of gravitational waves (GW) from astrophysical sources is one 

of the most outstanding problems in experimental gravitation today. Large 
laser interferometric gravitational wave detectors like the LIGO, VIRGO, LISA, 
TAMA 300, GEO 600 and AIGO are potentially opening a new window for the 
study of a vast and rich variety of nonlinear curvature phenomena. This network 
of gravitational wave detectors can confirm that GW exist and monitoring these 
O ' wave forms give important information on their amplitudes, frequencies and 

other important physical parameters. 

The detection of GW necessitates the substantial accumulation of Signal 
to Noise (S/N) over long observation periods. The data analysis for continu- 
ous GW like, for example, from rapidly spinning neutron stars is an important 
problem for ground based interferometric detectors that demands analytic, com- 
^ ' putational and experimental ingenuity. 

In recent works |2 El we have implemented the Fourier transform (FT) of 
the Doppler shifted GW signal from a pulsar with the Plane Wave Expansion in 
Spherical Harmonics (PWESH). Spherical-harmonic multipole expansions are 
used throughout theoretical physics. Indeed, they arise wherever one deals with 
fields, be they electromagnetic, gravitational, hydrodynamical and solid body 
etc. The expansion of a plane wave in spherical harmonics has a variety of appli- 
cations not only in quantum mechanics and electromagnetic theory I" , but also 
in many other fields. In linear theories, such as vacuum electromagnetic-wave 
theory, the problem is relatively simple: the field's multipole components evolve 
independently of each other; there is no coupling. However, in nonlinear theories 
like general relativity, it is more difficult. A number of researchers have used 
spherical-harmonic expansions for a variety of problems in general relativity, in- 
cluding problems where nonlinearity shows up as Kip Thorne |S]has pointed out 
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in detail. The relationships between scalar spherical harmonics incorporated in 
our present work and the various other harmonics like the pure spin vector har- 
monics related to the Regge- Wheeler harmonics have been discussed [S]. The 
basis states in the PWESH expansion form a complete set and facilitate such a 
study. It also turns out that the consequent analysis of the Fourier Transform 
(FT) of the GW signal from a pulsar has a very interesting and convenient de- 
velopment in terms of the resulting spherical Bessel, generalized hypergeometric 
function, the Gamma functions and the Legendre functions. Significant analytic 
and numerical studies have been carried out by many, for example, by Olver [H], 
R.C. Thorne [7|, van der Laan and Temme|H] and Cherry Hypergeometric 
functions also arise in the analysis of gravitationally radiating binary stars |1U|. 
The generalized hypergeometric functions in our analysis find their extensions 
to the H-functions, Meijer G-functions |11[I12| and the Heun |13[I14| functions. 

Both rotational and orbital motions of the Earth and spindown of the pul- 
sar can be considered in this analysis which happens to have a nice analytic 
representation for the GW signal in terms of the special functions above. The 
signal can then be studied as a function of a variety of different parameters 
associated with both the GW pulsar signal as well as the orbital and rotational 
parameters. Regardless of the immediate usefulness of such an analysis for GW 
data community which we do indeed look forward to the numerical analysis 
of this analytical expression for the signal offers a challenge for efficient and 
fast numerical and parallel computation. The plane wave expansion approach 
was also used by Bruce Allen and Adrian C. Ottewill in their study of 
the correlation of GW signals from ground-based GW detectors. They use the 
correlation to search for anisotropics from stochastic background in terms of 
the l,m multipole moments. Our PWESH formalism enables a similar study. 
Recent studies of the Cosmic Microwave Background Explorer have raised the 
interesting question of the study of very large multipole moments with angular 
momentum / and its projection m going up to very large values of I ~ 1000. 
Such problems warrant an intensive analytic study supplemented by numerical 
and parallel computation. 

However, the important spin down corrections were not included in our pre- 
vious works. A considerable amount of work has been done on the search for 
the continuous GW signals. Livas, Jones and Niebauer ^B] have investigated, 
in separate works, the time series of continuous GW signals that incorporated 
modulation of the motion of the detector over limited regions of parameter 
space. They have not considered pulsar spin down, and restricted their anal- 
ysis to small areas of the sky. In this work, we present a parametrized model 
clearly discussed by Brady et al P EI for the gravitational wave frequency and 
the phase measured at the ground based detector that includes the spindown 
parameters. 

Spin down corrections are important in the analysis of GW signals. Brady 
et al indicate that most of the S/N is accumulated during the final stages 

of spin down of pulsars. There could exist a class of pulsars which spin down 
mainly due to gravitational radiation reaction 18 . The frequency scales as 
/ oc r^^/* in these types of pulsars, where t is the age of the pulsar. Assuming 
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that the mean birth rate for such pulsars in our galaxy is , the nearest one 
should be a distance r = \/tb/t from earth, where i? ~ 10 kpc is the radius 
of the galaxy. This estimate has been provided by Brady et.al |T]- We use the 
expression for he provided by Thorne ^Hl and |20l 1211 1^ . 
The amplitude, he is given as, 

= — (1) 

where G is the gravitational constant, / is the sum of the frequency of 
rotation of the star and the precession frequency, / is the moment of inertia 
with respect to the rotation axis, e is the poloidal ellipticity of the star and r is 
the distance to the star. Hughes et. al, suggest that he < 10~^^. 

The frequency domain characteristic of the GW signal consists of two com- 
ponents with carrier frequencies /o and 2/o (Jaranowski et. al.l20') that are 
both amplitude and phase modulated. The amplitude modulation in their anal- 
ysis is determined by functions which split each of the two components into five 
lines that connect the GW and the earth rotation frequencies. They observe 
that the frequency modulation (FM) broadens the lines. They make estimates 
for a kilohertz GW frequency, the spin-down age r = 40 years, and an obser- 
vation time of 120 days for the maximum frequency shifts due to the neutron 
star spin-down. Earth's orbital motion and Earth's diurnal motion which are, 
respectively, of the order ^ 8, ~ 0.1, and ^ 10^^ Hz. In this paper, we have 
considered a single component GW frequency, /o and neglected the amplitude 
modulation. 

The paper is organized as follows. In Section 2 we outline the FT of the 
GW signal from the pulsars which is relevant to the detection of continuous 
gravitational waves. We focus attention on spin down effects which indicate 
significant frequency evolution over periods of several weeks of observation. In 
this section, we use the expression for a parametrized model of the expected 
gravitational waveform, including modulating effects due to detector motion. 
We derive the analytic form of the FT with spin down corrections in the following 
section. The expressions for the spin down moments are computed using the 
familiar technique of differentiation with respect to a parameter. We present 
our conclusions in section 4. 



2 Frequency and phase evolution with Spin Down 

Manchester (1992) and Kulkarni (1992) have suggested that pulsars lose ro- 
tational energy by electromagnetic braking, particle emission and emission of 
GW |23[ I24| . Therefore, the rotational frequency is not completely stable, and 
varies over a timescale t which is of order the age of the pulsar. Younger pulsars 
with periods of tens of milliseconds have the largest spin down rates. Current 
observations suggest that spindown is primarily due to electromagnetic brak- 
ing. Brady et. al. have suggested a sufficiently general model of the frequency 
evolution to cover all possibilities in the study of pulsar GW signal detection. 
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For observing times tabs much less than t, the frequency drift is small. We use 
a slightly modified form of the parametrization given by Brady and Creighton 
(2000)^7] in that it differs by an factor of 1/2, and the summation index, k 
starts from zero. Our expression is the Doppler modulated version of the earlier 
power series of the frequency evolution given by pp. The parametrized model 
for frequency and phase evolution with spin down corrections is given in the 
following form. 



/(i;A) = 



^(i; A) = tt/o 




(2) 



(3) 



where r(t) denotes the detector position, h — (sin cos 0, sin sin i^, cos 0) 
is a unit vector in the direction of the source. The angles 6, (j) are associated 
with the pulsar source; v{t) = ^^^^ is the velocity of the detector and ft are 
the spin down parameters. Here \fk\ < T^in ^^nd Tmin is the spin down age in 



years. Here r = cif /dt ^^^'^ years of the pulsars. For the extreme case of 
the gravitational-wave frequency of 10'^ Hz, the spin-down age is r = 40 years 
for pulsars in our galaxy and > 10^ years for millisecond pulsars. Our analysis 
considers the equality situation, that is, |/fc| — t~^„. The phase (j) depends 
on the frequency /o, k spin-down parameters, and on the angles (f>, 9 and the 
CO- latitude a of the detector. 

In the earlier analysis of Valluri et.al. [3], the spin down corrections were 
neglected. In their analysis of the FM of a monochromatic plane wave, they 
characterize the motion of the Earth (and detector) by: (a) assuming the orbit 
of the Earth to be circular, (b) neglecting the effect of the Moon and the 
perturbation effects of Jupiter on the Earth's orbit. In this work, in addition to 
incorporating spin down effects, we also include the approximate correction for 
the eccentricity of the orbital motion of the Earth. 

With these assumptions, (a) and (b), and ignoring spin down corrections 
we arrive at the following expression of the FT for the GW pulsar signal P], 
rewritten more concisely as. 



oc oo Z 

h{f) = S'„im(wo, ^orb, Tev, n, I, m, A, i?, k,a,6,(j))= X! X! X! V'oV'iV'2V'3V'4 

n= — ooi=Om= — i 

(4) 

where 



i}:o{n, I, m, a, 6, cj)) = Am^Yi„,{e, 0)iV;„jF;"(cos a) 



(5) 
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(6) 



V'4(fc, ^, n, Uorb, t^r) =1 F3 N + 1; Z + -, 



3 ^ + i?orb + 2 ^ — Sorb + 2 — fc^ 



2' 2 ' 2 '16 



(9) 



The angle a is the co-latitude detector angle and angles 0, (j) are associated 
with the pulsar source. Here loq — 27r/o, ijJorb — 5^ {Torb ~ 365 days, Tsr — 1 

day), Borb - 2 + f + ^^), A: = 47r/oflE sin(a) ^j^g ^^^.^g 

Earth, c is the velocity of light) and A is the sun-earth distance. 



3 Calculation of the FT including Spin Down 
Corrections 



Considering the spin down corrections in Equation (2), we have for the summa- 
tion term in (j){t; A) in Equation (3), 



fe-i-i 



Upon using the definition, \fk\ = Tmin have for this summation, 



fc+i 



This power series represents a logarithmic function of the form — ln(l — x). It 
is interesting to note that this is the n — 1 case of the multi- valued polylogarithm 
function Ln{x) and the phase evolution, apart from the factor iTrfo is one of the 
versions of the multi- valued Lambert- W function 12^. Hence, we have, 



In 1 



cj){t) =exp[z7r/o {(t + ^) - T„„„ln(l - ^)} 



e exp 



"«T™„j7r/oln(l - -7—^) 
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1 - 



(10) 



Using the binomial expansion on equation (10), we obtain, 



1 + infoT„ 



{iT^ foTrmn){™ foTmin + 1) M + 



We define the Spindown Moment Integrals as follows: 



V c 



dt 



(11) 
(12) 



for k — 0,1,2,3,...., where fc = denotes the absence of spin down corrections 
and recovers the formula for the FT given in Equations (4-9). 

Thus, a generic spin down moment integral could be written from the deriva- 
tive of a generic integral. 



^generic 

Hence, we have the /c— th derivative as 
T 

^ ^generic 



(i7r/o(t+£j2:)-i27r/t)^^ 



(13) 



for k =0,1,2,3,.. 













r ■ n\ 


f 


in -\ 




Jo 







This k derivative gives the Corresponding Spin down moment integral. Thus, 
Igeneric has bccu analytically evaluated. Its derivatives with respect to /o can 
be done by symbolic packages like Maple, Mathematica etc. It should be noted 
that the FT of the GW signal using the later parametrization given by Brady 
and Creighton (2000) can also be analytically incorporated in a similar manner 
of derivation given above. 

It should be noted that previously these spin down corrections were incor- 
porated by splitting up a time interval, say Tq into M equal parts each of 
interval At {Tq = MAt) so that the signal is monochromatic in each interval or 
"window" [22]. This is the idea behind "stacking" where FTs are incoherently 
combined by adding their power spectra. "Tracking" as suggested by Brady 
et. al (2000), Papa and Schutz, (cited in review article gr-qc/9802020) is the 
attempt to track weak lines in individual FTs in successive data sets to identify 
persistent signals. We express our final result as a coherent sum over any given 
time interval as a closed form solution. 



In our original FT, we have not included the effects of orbital eccentricity 
which arise due to the orbital motion of Earth. Denoting the orbital eccentricity 
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by 60 (= 0.017), we find that by replacing the usual Sun-Earth distance, A by 
the following modified A, is a reasonable estimate of the correction. 

A ^ A^l-el (15) 

The correction in the phase of the GW signal due to the effects of orbital 
eccentricity can also be accurately estimated as has already been pointed out 
by Jaranowski et al. (1998). 

4 Conclusions 

Recently, new mechanisms e.g. r-mode instability of spinning neutron stars and 
temperature asymmetry in the interior of the neutron star with misaligned spin 
axis have been discussed in the literature [21 122 • The continuous G W signal 
may consist of frequencies which are multiple of some basic frequencies. Brady 
and Creighton have given estimates of the characteristic strain of gravitational 
waves from an active r-mode instability in a newborn neutron star suggesting 
that these sources will be detectable by the enhanced interferometers in LIGO 
out to distances 8 Mpc; the rate of supernovae is 0.6 per year within this 
distance plj. [T7j. 

We have presented in this paper the rudiments of a simple analysis for spin 
down for sources of continuous gravitational waves. For the more computationally- 
intensive search over all sky positions and spin down parameters, it is impor- 
tant to be able to calculate the smallest number of independent parameter 
values which must be sampled in order to cover the entire space of signals. The 
PWESH has the potential to improve the numerical accuracy and convergence 
of analytic FT's and spin down corrections associated with the GW signal. 

The study of templates in time-domain has been made by many research 
workers(Schutz 1991; Krolak 1997; Brady et al. 1998; Brady, Creighton 2000; 
Jaranowski et al. 1998; Jaranowski, Krolak 1999, 2000 26 ). However, the 
analysis in the frequency domain has the advantage of incorporating the spectral 
noise density of the interferometer. The data output at the interferometer is 
available in discrete form and the question arises if the analytical FT is in fact 
a very convenient tool. However, the FFT has a resolution limited to I/Tq as 
pointed by Srivastava and Sahay [221. hoped that the analytical FT will 

provide useful insights. 

We intend to further estimate from our spin down analysis the number of 
templates required for matched filtering analysis. We intend to make use of 
our approach developed for all frequencies to evaluate the Fitting Factor (FF) 
as given by Apostolatos |2Z1 and considered for low frequencies by Srivastava 
and Sahay (2002). The time-delay type effects |2S| due to general relativity 
and other corrections for the ephemeris have been treated in detail by Cut- 
ler in his Ligo Algorithmic Library (LAL) Barycenter Package (available at: 
|http: / /w ww.lsc-group.phys.uwm.edu/ ) and will also be taken into account. 

A new expansion was explored by MacPhie et.al. ^ both in the case of scalar 
and vector harmonics. This can be extended to tensor harmonics as well, and 
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work in that direction will be carried out. This would be of possible relevance 
to study the multipole formalisms for gravitational radiation |S] . 
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